Can the dark matter halo be a coUisionless ensemble of axion stars? 
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fa > lO^GeV |14l423l |. Cosmological observations im- 
pose constraints on nia to be bigger than 10~^eV and 
fa < lO^^GeV f23 - [27t . Therefore, cosmological and as- 
trophysical observations put a range on to^ to be 



< 77la < lO^^cV. 
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If dark matter is mainly composed of axions, the density distribution can be non uniform dis- 
tributed but clumpy instead. By solving the Einstein-Klein Gordon system of a scalar field with 
a potential energy density of an axion like particle, we obtain the maximum mass of the self grav- 
itating system made of axions called axion stars. The collision of axion stars with neutron stars 
may release the energy of axions due to the conversion of axions into photons in the presence of the 
neutron star magnetic field. We estimate the energy release and shown that it may exceed the solar 
luminosity per collision but should be much less than previous estimates 0, 0I- Future data from 
femtolensing should strongly constrain this scenario. 
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I. INTRODUCTION 

Nowadays one of the most interesting problems in 
physics is to reveal the nature of the dark matter (DM). 
The existence of DM is supported by large amount of 
astrophysical observations, in both galactic and cosmo- 
logical scales d, Q ■ Since those effects are purely grav- 
itational, the standard candidate is a massive particle 
with null or very weak interaction with the rest of the 
particles of the standard model of particles (SM). Two 
of the leading candidates for dark matter are the ax- 
ions and weakly-interacting massive particles (WIMPs) 
as neutralinos Q- 

The axion is the pseudo Nambu-Goldstone boson gen- 
erated in the spontaneous breaking of the U{1) Pcccci- 
Quinn (PQ) global symmetry PQ symmetry was 

introduced to explain the smallness of the strong CP vio- 
lation in QCD. Initially, this new field was assumed to be 
"invisible" as very weakly interacting but in 1983, it has 
been shown that axion couples to two photons through 
PrimakofF effects @ and could be detected through its 
conversion into a photon in presence of a strong mag- 
netic field [13. 

Most of the direct detection methods for the axion are 
based on the conversion of axions into photons [ll| - [l3| . 
But it is also possible to constrain the axion mass and 
coupling through astrophysics observation. As usually, 
stars have strong magnetic fields. Sun or any other as- 
trophysical object as dwarf stars, supernovae, etc. could 
produced axions in their core during cooling processes 
and to transform them into photons through PrimakofF 
effects. The search for photons produced in stars through 
axion conversion can be used to put an upper limit on 
axion mass nia and coupling fa where fa defines the 
scale of U{l)pQ symmetry breaking: iria < 10~^eV and 
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So, an axion mass within this r ang e could be a nice can- 
didate for Cold Dark Matter [28l [29| as it could form 
Bose-Einstein condensate [30l |. 

The rate of detection needs a prescription on the dark 
matter distribution in the galactic halo, and in particu- 
lar, the local dark matter density. There are arguments 
in favor of a non uniform distribution of dark matter in 
the galaxy [sil - ls^ . For instance, early fluctuations in the 
dark matter may go nonlinear long before photon decou- 
pling, producing density-enhanced dark matter clumps. 
The evolution of the axion field at the QCD transition 
epoch may produce gravitationally bound miniclusters 
of axions [H, |3^. Such minicluster, due to collisional 
2a 2a process, may relax to a selfgraviting system 
[37I . [ssj . Those self-gravitating systems made of axions 
are simply called axion stars (AS). Previously, some au- 
thors explored the possibility that such axion stars may 
have some observable consequences. For instance, the 
collisions between axions stars and neutron stars may 
indu ce g amma ray burst [39| . ultra- high energy cosmic 
rays [40|, or induce emission of X-ray by neutron stars 
[H . All those observables have in common the possibility 
that the axion may oscillate into a photon. 

In this paper, we study the self-gravitating system 
made of axions with a realistic potential for the scalar 
field as in [4l| contrary to [11, [33, \^ . In section II, we 
obtain the typical Mass-density curve varying the cen- 
tral density. This curve exhibits a maximum mass which 
depends on the value of the axion mass and the central 
density. Section III is dedicated to update signatures of 
the AS in the Galactic Halo due to AS collisions with neu- 
tron stars (NS). The huge magnetic field in NS, B - 10^° 
Gauss, may convert to photons all axions involved in the 
collisions. In this section, an estimate of the total radi- 
ated energy is computed. Finally, we present our results 



and compared to previous estimates in Section IV. 



II. AXION STARS 

A. On the formation of axion stars 

The formation of axion miniclusters were studied in 
[35I [37I Issj . They did a numerical study of the evolu- 
tion of a spherically symmetric axion fluctuation in an 
expanding universe at the QCD epoch, where the mass 
of the axion effectively switches on. The axion potential 
considered was 



^(0) = r>^lfa 
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where (f> is the axion field, fa the energy scale of decou- 
pling, and rua the axion mass. The study was done for 
a wide range of initial conditions and a common feature 
of the final density distribution is that it developed a 
sharp peak in the center, i.e. due to the attractive self- 
interaction, some non-linear effects induce axion pertur- 
bations that lead to very dense axion miniclusters. If such 
density peaks are high enough, gravitational collapse and 
subsequent virialization can occur. This relaxation is due 
to 2a — > 2a scattering and the associated relaxation time 
is smaller than the age of the universe if the energy den- 
sity of the minicluster is p > 10^°eV^ [13]. An improved 
numerical resolution of the axion mini-clusters revealed 
that around ~ 10% of the axion density will be in mini- 
clusters with densities that can virialize to coherent axion 
fields on a self-gravitational well, i.e. axions stars [4^ . 

Other mechanism that may lead the formation of axion 
miniclusters and dark matter miniclusters was proposed 
in (36j . There, at the time where the phase transition 
from a quark-gluon plasma to a hadron gas take place, 
the spectrum of density perturbations develops peaks and 
dips produced by the growth of hadronic bubbles. Af- 
ter wards, kinematically decoupled cold dark matter falls 
into the gravitational potential wells provided from those 
peaks leading the formation of dark matter clumps with 
masses < lO'^^M^. 



8e-21 



1 6e-21 - 



4e-21 



O 



2e-21 - 




5e-15 



3e-15 



2e-15 







J , , L 



2e+ll 4e+ll 6e+ll 3e+10 6e+10 9e+10 

p(0)[Kg/m'] p(0)[Kg/m^] 



FIG. 1: Mass of axion star. Left Panel: ma — 10 ^eV. Right 
Panel: rUa = 10"^eV. 



where the source of Einstein equations, the energy mo- 
mentum tensor (T^,^), is the average over the ground 
state of a real, quantized scalar field ^{r,t) with poten- 
tial given by the energy potential ([2]). We work in a 
spherically symmetric metric ds'^ = B{r)dt^ — A{r)dr'^ — 
r2(sin^ edcjp'+de) , and to assume a harmonic dependence 



of the field $(r,i) 
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to the following system of equations: 
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B. Properties of the axion stars 

The previous picture on the formation of axion stars 
did not study the final self-gravitating system made of 
coherent axion fields. In order to do so, we have solved 
the Einstein-Klein-Gordon equation in the semiclassical 
limit 
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where we have defined 



fa 



B = 



B 



fa V 47r ' rria 

. In the present work we extend the results of [4lj by 
solving the same system but for a wider range of masses 
of the axion field and many initial values of the axion field 
at the center cr(x ~ 0). Since the axion mass is restricted 
to be 10~®eV < iria < 10~'^eV, we show in Fig. [T] the 
resulting masses of the AS as a function of the value of 
the central density p(0) = to^(t(0)2/2. The mass-density 
plot shows that there exist a maximum mass for a par- 
ticular value of the central density. The density profile of 
those maximum mass configurations are shown in Fig. [2] 
and the numerical values arc reported in Table ID Those 



, as it was done in 
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FIG. 2: Density profile. Top Panel: = 10"^eV. Bottom 
Panel: rUa = 10"^eV. 
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FIG. 3: Number of axion stars with a mass of 10 AIq for 
different dark matter profiles. 



TABLE I: Masses and typical radio Rgg for the maximum 
mass configurations shown in Fig. [T] where Rgg is the radius 
of the AS where 99% of the mass is contained 



Axion 


mass 


(eV) 


p(0) (Kg/m=^) 


Mass (Mq) 


Rgg (meters) 


rria 


= 10" 


-5 


2.1 X 10^° 


5.0 X 10"^^ 


119.40 


ma 


= 10" 


-3 


6.8 X 10^° 


7.4 X 10"^^ 


0.89 



configurations arc our benching mark in order to study 
possible astrophysical signatures. 



III. POSSIBLE AXION STAR SIGNATURES IN 
THE GALACTIC HALO 

The dark matter halo of galaxies can be considered as 
a coUisionless ensemble of mini-MACHOs, as far as this 
hypothesis is not in contradiction with the limits imposed 
by microlensing or gravothermal instability. In particu- 
lar, scalar field mini-MACHOS has been considered as 
a possible realization if their masses are below IO^'^A/q 
[43| to evade the microlensing limits. In turns out from 
our analysis on the maximum mass for an axion star, for 
nia = 10~^eV that axions stars may play the role of such 
mini-MACHOS. In addition to microlensing, femtolens- 
ing could be a useful tool in order to look for axion stars 
[42l | . Recently, new limits from femtolensing of gamma 
ray burst provides new evidence that primordial black 
holes in the mass range 10^^^ — 10^^^ Mq do not consti- 
tute a major fraction of the dark matter [i^l and once 
again, the maximum mass of the axion star are consis- 
tent with those limits. In the case for an axion mass 
of nia = 10~^eV, the maximum AS mass is far to be 
excluded from these current limits from femtolensing of 
GRBs since M^ax ^ lO^^^M© but it is encouraging that 
in the future, femtolensing measurement may confirm or 
deny the possibility that the dark matter halo may be 
made of an ensemble of axion stars. In addition to AS, 
other types of dark matter clumps may play the role of 



mini-MACHOS. In particular if dark matter is a non self- 
annihilating fermion, for some values of its mass it may 
fulfill the requirements to be a mini-MACHO [45|. For 
a self-annihilating fermion as the ncutralino, neutralino 
stars j46l - l48| . have been proposed. Actually, the max- 
imum mass for a neutralino star is ~ 10"'' making 
them suitable dark matter mini-MACHOS candidates, 
although their stability is questionable (ioj . 

The clumping of dark matter may change the expected 
number of events for direct dark matter searches since the 
local density may be largely affected at the vicinity of the 
Earth. Nevertheless, the clumpy dark matter may offer 
new ways for indirect detection. In the case of the axion, 
early studies of luminous axion clusters where done in 
the 80 's 5C|. In this work we reconsider the proposal of 
Iwazaki |39j where axion stars interact with the strong 



magnetic field of a neutron star producing photons, but 
now we consider the mass and radius of the AS obtained 
as general relativity indicates, by solving eqs. [5] 

We can estimate the number of AS in a galactic halo 
using dark matter density profiles. In Fig. [3] it is shown 
the number of axion star with a mass of 10~^^Mq for four 
different DM density profiles: the Navarro Frenk White 
profile 51 1, the Moore [s^], the Kravtsov [H^, and the 
modified isothermal profile [s^ ]. Notice the high num- 
ber of axion stars at the center of the galaxy. Further- 
more, it is expected that a galaxy has around ^ .1% of 
its total stars will finish their life as neutron stars. The 
Milky- Way, for instance, could have around ~ 10^ neu- 
tron stars. The high number of axion stars (see Fig. [3]) 
and the expected number of neutron stars implies a non- 
zero probability that AS collide with a neutron star and 
dissipate their energy under the effect of the neutron star 
strong magnetic field. 

In order to estimate the amount of energy radiated by 
this process, we closely follow [s^. Starting with the 
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FIG. 4: Energy radiated by the collisions of axion stars with 
neutron stars in the dark matter halo of galaxies 

Lagrangian 

L = Ud^ad^a-m^c?) - ~F^,F^- - -F^.F^^ , (6) 

here c is a coupling constant of order unity, a is the axion 
field and F^^'^ is the electromagnetic stress tensor and F^'" 
its dual Ell . And expressing this Lagrangian ^ in 
terms of the electric and magnetic fields E^B. 

^aii = -r^F ■ B (7) 

Ja 

we can now derive a "modified" Gauss law, namely: 

OE^^d-iaB). (8) 

Ja 

Thus, the axion field a may induce an electric field 
Ea = —caaB / fa- If the neutron star has an electric 
conductivity a, an electric current Jm ~ <jEa is pro- 
duced. The dissipation in the magnetized conducting 
media, with average u electric conductivity is evaluated 
using the Ohm's law W = aE^cfix, which, with help 
of the modified Gauss Law, can be directly related with 
the axion field and the neutron star magnetic field: 

W = 47r^^^ r a{r fr^dr. (9) 

Ja Jo 

Taking the two profiles for the axion field shown in Fig. 
m we perform the integration of eq. ([9]) numerically. The 
radiated energy is shown in Table HIl Finally, we estimate 
the rate of collisions, such as we estimate the luminosity 
radiated every second. The number of collisions per pc^ 
per second will be 

Rc = riAsir) X pNsir) x S xv, (10) 

where uas is the number of AS per pc^ as a function of 
the distance to the galactic center (Fig. [T]), pjys is the 



TABLE II: Energy radiated for an AS with the maximum 
mass configurations shown in Fig. [2] for a NS magnetic field 
B = 10*G and an electric conductivity a = lO'^'^s'^. [13] • 



Axion 


mass 


(eV) 


fa (GeV) 


Mass (Mq) 


W (Lq) 


rria 




-5 


6 X 10^ 


5.0 X 10"^^ 


1.6 X 10^° 


rria 


= 10" 


-3 


6 X 10" 


7.4 X 10^^^ 


2.4 X 10* 



probability to find a Neutron Star at that point, and we 
will assume this distribution to be [56j 

pjvs(r-) = ^r"-VA"e-'-/\ (11) 

where a ~ 1.7 and A(kpc) = 5.21 as given in ref. [56j . 
For simplicity, the z— dependence is neglected. This as- 
sumption implies that NS distribution is assumed to be 
spherically symmetric. A is the normalization constant 
defined such that 

10^NS = J ApNs{r)r^sm^9drded(l}. (12) 

Finally, S is the cross section which is simply computed 
as the effective area of the AS, i.e. S = ttRIq. The 
velocity v of AS is assumed to be around V = 2xl0^cm/s. 

The energy radiated per second is simply L = W x Rc 
which is shown in Fig. 2] for the case of the maximum 
mass of the axion star M = 1O~^^M0. 



IV. CONCLUSIONS 

As we can see from Table HIl the energy dissipated by 
one collision of an AS with a NS is many orders of mag- 
nitude bigger than the solar luminosity. But our results 
for the total dissipated energy per second in NS-AS col- 
lisions are around 10 orders of magnitude less than pre- 
vious estimates P, [11 • As it can be seen from Fig [TJ this 
difference comes from our computation of the maximal 
AS mass obtained through resolution of EKG equations. 
The obtained maximal AS mass is much smaller than the 
corresponding parameter used in ref. [H [H ■ It is impor- 
tant to notice that the maximal AS mass used in these 
references are now excluded by last observational results 
from Femtolensiong of Gamma Ray Bursts (i^ . 

Furthermore, improvement in femptolensing will help 
us to constrain the axion mass by constraining the max- 
imum mass of AS. Observation of a sudden release of 
energy in Neutron stars may indicate the interaction of 
an AS with its strong magnetic field. If some of those 
signatures are observed, they can be considered as indi- 
cations that the dark matter halo is made of an ensemble 
of axion stars. In our estimations, we have considered the 
largest mass of the AS and a hypothetical electric con- 
ductivity a [s^!, nevertheless, we hope our results will 
encourage further studies in the possibility of detecting 
the axion by its possible clustering. 
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